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Abstract. We give two derivations of magnetic flux quantization in a 
superconducting ring in the shape of a Mobius band, one using direct 
study of the Schrodinger equation and the other using the holonomy 
of flat (7(l)-gauge bundles. Both methods show that the magnetic flux 
must be quantized in integral or half-integral multiples of 4>o = hc/(2e). 
Half-integral quantization shows up in "nodal states" whose wavefunc- 
tion vanishes along the center of the ring, for which there is now some 
experimental evidence. 



One of the best-known macroscopically observable quantum effects is the 
quantization of magnetic flux through a superconducting ring M in units 
of $ = || [2]. The question we treat here is this: What happens to this 
condition if the superconducting ring M is in the shape of a Mobius band 
(Figure [T]) rather than an annulus? Does the condition remain the same, or 
should it be modified? This provides a good test case for the application of 
topology in physics. It should be possible to check the result experimentally 
since crystals of potentially superconducting materials such as NbSe3 have 
recently been produced with a Mobius band shape O [10] . Investigation of 
a different topology, the "figure eight," was carried out by Vodolazov and 
Peeters [TT] . 

Several possible derivations for flux quantization have been given. We 
shall analyze what happens to them under this "exotic" topology, under the 
assumption that the superconductor is extremely thin and can be treated 
as if it were 2-dimensional. We can take coordinates in the superconductor 
to be x and y, with x (going around the ring M) going from to y going 
from — w to w, and the point (0, y) identified with (£,—y). The simplest 
treatment (following [7], based roughly on the classic paper [I]) starts with 
the Schrodinger equation for a Cooper pair of charge — 2e in an electromag- 
netic field with circular symmetry. Recall that the vector potential A is 
not necessarily globally defined on the ring, though we can take it to be 
single- valued on [0,^] x [— w,w], possibly with different values at (0,y) and 
at (£, —y). By the Meissner effect, B = VxA = 0in the superconductor 
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Figure 1. A Mobius band with its coordinate system 

M, which implies that A is locally of the form V</>. Again we can take (f> to 
be single- valued on [0,£] x [— w, w], possibly with different values at (0, y) 
and at (£, — y). The (time- independent) Schrodinger equation becomes 

(0.1) — (-iHV+ — A) 2 ^ + FV = E*I>, 

2m c 

where is the wavefunction of a Cooper pair, treated a single boson with 
mass m. Let be a solution of (jO.ip with A = 0. Then if ip = ijjQe -10 "^ for 
a constant a, we have 

(-ihV + ^A)-0 = (-iW^o)e - ^ - toVoe" iQ0 A + ^oe" iQ0 A 
= -ihe-^Vipo, 

provided that a = ||. Thus ^ will be a solution of (jO.ip for this value of a. 
Since the wavefunction must be well-defined globally, in a ground state where 
ipo is everywhere non-zero on the interior of M, (p(l, —y) — 1^(0, y) must have 
a constant value of 27r 2 ^ c = = n<E>o> for some n G Z, independent of y. In 
particular, if C is the closed curve through the middle of the superconductor, 
corresponding to y = in our coordinate system (see solid black curve in 
Figured]), then <f c A • dr = n&o, which is the flux quantization condition. 
Note however that if we take the curve C given by y = const, with the 
constant non-zero, then when x runs from to £, we end up on the opposite 
side of C from where we started, and so one has to let x run all the way out 
to 21 to get back to the starting point (see dashed black curve in Figured]). 
Thus for C', we have the modified quantization condition <f c „ A-dr = 2n&Q, 
with an extra factor of 2. 
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However, another phenomenon, proposed on different theoretical grounds 
(calculations with the Ginzburg-Landau and Bogoliubov-de Gennes models) 
in [U [5], and supported by experimental evidence in [6], is also possible. 
Namely, one can have a "nodal" state supported near the dashed curve C, 
with i/>o = on the circle C given by y = 0. Then it suffices to for <f)(£, y) — 
4>(0,y) to be a half-integer multiple of <I>o- This will still give a globally 
single-valued wavefunction since the point (£, y) is identified with </>(0, —y) 
and <f>(£, —y) — </>(0, —y) is again a half-integer. Thus the flux appears to 
be quantized in half-integral units of as was observed in [6]. A slightly 
different analysis by Mila et al. [8] showed that for a Mobius "ladder" 
geometry, quantization should be in multiples of $o " as long as coherent 
motion between the chains is possible," and in multiples of Q$/2 when there 
is no coherent motion between the chains (so that particles effectively travel 
along the curve C). 

A more sophisticated approach follows the ideas of [7]. Following Dirac 
[3j and Weyl, we view the magnetic field as a J7(l)-gauge field. More pre- 
cisely, the vector potential A is a connection on a £/(l)-bundle and the field 
strength (curvature 2-form) is the magnetic field. The phase of the wave 
function ip is a section of this £7(l)-bundle [3]. Because of the Meissner ef- 
fect, the curvature of the bundle vanishes on M, i.e., the bundle is flat. Since 
H 2 (M,7j) = 0, the bundle is also topologically trivial, and the connection 
differs from the usual connection (corresponding to the case A = 0) by a 
constant 1-form, which we can identify with a single real number, which is 
the value of the holonomy or flux $ = <f c A-dr. Since the wavefunction 
must be single-valued on M, and since as we saw above, the change in the 
wavefunction as we go around the loop C is 

e ft C , 

we obtain the flux quantization condition = n&Q. 

The "nodal state" case can be treated similarly, except that we replace 
M by the complement of C since we are assuming the Cooper pairs are 
localized away from the center of the Mobius band. (If the wave function 
vanishes on C, then its phase there is not well-defined, so the bundle is not 
defined on C.) Since the inclusion (M \ C) ^ M sends a generator of the 
fundamental group of M \ C to twice a generator of the fundamental group 
of M, we have in effect twice as many possible flat bundles. In other words, 
we only require § c , A • dr to be integral, which means that flux around the 
superconducting loop satisfies the flux quantization condition <3? = n&o with 
n a half- integer. 

The second approach also explains the answer to another question: the 
flux $ = <f c A • dr only depends on the homology class of the loop C in 
Hi(Al,Z) = 7Ti(M) = Z, since this is a basic fact about holonomy of flat 
connections. In particular, since C is homologous to 2C, the flux around 
C is quantized in twice the units of the flux around C . From the first point 
of view, this is a bit harder to see, since if C\ and C2 are homologous loops, 



-1 



JONATHAN ROSENBERG AND YEHOSHUA DAN AGASSI 



the usual argument (when M is an annulus) would be to take the region D 
with boundary C\ — C2 (i.e., C\ U C2, but with reversed orientation on C2) 
and to use Stokes' Theorem to argue that 

(0.2) 0=//(VxA)-n=i A • dr - <b A -dr. 

J JD JCi Jc 2 

This runs into difficulties since M is not orientable (so that the normal n 
is not well defined), and thus Stokes' Theorem doesn't seem to apply. But 
one can rectify things by cutting M open along a curve whose complement 
is orientable. 

An issue not settled by our analysis is what happens with a "thick" Mobius 
band (thickness being measured compared to the penetration depth), where 
2-dimensionality is no longer a good assumption. In this case, since there 
may not be nodal states in the sense described above, one might expect flux 
quantization only in integral multiples of <J?o- It would be interesting to test 
for this experimentally, as the results of [6] are only for a thin film. 
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